Non-Markovian Relaxation of a Three-Level System: Quantum Trajectory Approach 
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The non-Markovian dynamics of a three-level quantum system coupled to a bosonic environment 
is a difficult problem due to the lack of an exact dynamic equation such as a master equation. 
We present for the first time an exact quantum trajectory approach to a dissipative three-level 
model. We have established a convolutionless stochastic Schrodinger equation called time-local 
quantum state diffusion (QSD) equation without any approximations, in particular, without Markov 
approximation. Our exact time-local QSD equation opens a new avenue for exploring quantum 
dynamics for a higher dimensional quantum system coupled to a non-Markovian environment. 
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Every small quantum system, such as a two-level atom 
(qubit), a three- level atom (qutrit) or a cavity mode (sim- 
ple harmonic oscillator), should be regarded as an open 
quantum system due to the inevitable interaction with 
its environment. The effect from the environment may 
bring about some grave problems, such as quantum de- 
coherence and disentanglement, to the dynamics of the 
open system pj. Much has been studied in the Markov 
limit where the dynamics of open quantum system is typ- 
ically described by a standard Lindblad master equation 
which is equivalent to a quantum state diffusion (QSD) 
equation for a pure state [2J47J (setting h = 1): 

j t i> t = -iH^ t + A t (L)(z t + (L^ t )^ t ~^A t (L^L)^ u (1) 

where the notation At (A) = A — (A) t for any operator 
A and (A) t = (ipt\A\ipt) denotes the quantum average. 
L is the system operator coupled to the environment, of- 
ten called Lindblad operator. The dynamics of the open 
system is driven by the classical stochastic process Zt- 
The reduced density matrix for the system of interest 
can be recovered by averaging quantum trajectories gen- 
erated by the QSD equation ©: = M[\ip t (z)){ip t (z)\]- 
Here M[-] denotes the ensemble average over the classi- 
cal noise. Besides many appealing features exhibited by 
pure state quantum trajectories, Eq. ([1} provides a very 
efficient numerical tool in solving quantum dynamics of 
Markov open systems. 

Non-Markovian environments have become increas- 
ingly important in recent times due to their relevance 
in explaining new experimental advances in high-Q mi- 
crowave cavities, photonic crystals and atom laser in 
BEC [§-14 1 . It is also evident from the recent progress 



in quantum information processing, that environmental 
memory may be utilized to generate or modulate entan- 
glement evolution of an open quantum system [TBI [l6| . 
Clearly, an approach that is capable of dealing with non- 
Markovian quantum system is highly desirable. A non- 
Markovian QSD equation for a general non-Markovian 
open system by Diosi, Gisin and Strunz has provided 



a powerful tool in dealing with the exact dynamics 
of quantum open systems irrespective of the coupling 
strength, the correlation time and the spectral density 
of the bosonic environments 17H24|. Despite extensive 



research, the explicit non-Markovian QSD equations only 
exist for a single two-level system, the quantum Brown- 
ian motion model and optical cavities due to the intricate 
functional derivative appearing in the fundamental QSD 



equation [171 lla 125- 28] . Clearly, the power of the non- 



Markovian quantum trajectory method cannot be read- 
ily unleashed unless the stochastic QSD equation can be 
cast into a convolutionless form in which the numerical 
simulations and analytical applications can be easily im- 
plemented [iistj. 

In this Letter, we present, for the very first time, an 
exact time-local QSD equation for the three-level dissipa- 
tive dynamics in the framework of non-Markovian quan- 
tum trajectory at zero temperature. Our treatments of 
the three-level model have opened up a new route of ex- 
ploring dynamics of higher dimensional quantum open 
systems. Notably, our method of establishing a time-local 
non-Markovian QSD equation can be applied to multiple 
qubit systems to deal with the time behavior of them 
such as the estimation of non-Markovian entanglement 
evolution and fidelity [31. 32 [. 

We consider a Caldeira-Leggett-like model [2§| involv- 
ing a system with Hamiltonian H describing a three-level 
atom, coupled linearly to a general boson environment 
consisting of a set of harmonic oscillators a^a^ (e.g., 
cavity modes). The total Hamiltonian for the system of 
interest plus environment can be written as: 



fftot = H + ^(^LW + gicLal) + 



^w k a k a k , (2) 

k 



where the three-level system Hamiltonian H = loJ z = 
w(|2)(2| - |0)(0|), the Lindblad operator L = J_ = 
\/2(|0)(l| + |1)(2|) and = J + . Note that u is the 
spacing of the two neighboring energy levels of the three- 
level system, are the coupling constants between the 
three-level system and the environmental modes. It is 



2 



known 17j that the linear QSD equation for the three- 
level system can be formally written as: 



dt 



ipt{z) = -iujJ z 4> t (z) + j-z t ip t {z) 
dsa(t, s) 

0Z S 



J. 



(3) 



where a(t, s) is the bath correlation function and z t 



is a complex Gaussian process satisfy- 



ing M[z t ] = M[z t z s ] = and M[z* t z s \ = a(t,s). When 
a(t,s) = 5{t — s), it reduces to the memoryless Markov 
case. The density matrix of the system is recovered from 
the ensemble average of many realizations of quantum 
trajectories: 



Pt = M[\Mz))(Mz)\] 



(4) 

One advantage of using quantum trajectory is that, 
while deriving non-Markovian master equation is known 
to be a daunting task, the QSD equation governing the 
non-Markovian open system can be read off directly once 
the Hamiltonian of the total system is given by ([2} . How- 
ever, the formal non-Markovian QSD equation © or 
its nonlinear version (see below) cannot be easily im- 
plemented as an analytical or numerical tool unless the 
term containing functional derivative can be written as 
j^-ipt(z) = 0(t,s,z)ip t (z), where 0(t,s,z) is an oper- 
ator acting on the system Hilbert space satisfying the 
initial condition 0(s, s, z) = J_. The equation of motion 
governing the O-operator (consistency condition [3]) is 
given by 



d t O(t, s, z) = [—iujj z + J-z t — J + 0(t, z),0(t, s, z)] 
SO(t,z) 



5z s 



(5) 



where 0(t,z) = J Q dsa(t, s)0(t, s, z). For the model 
given by ([2]), it can be shown that this nonlinear operator 
equation can be solved by the following ansatz solution: 

O = f(t,s)J-+g(t,s)J z J-+i f p(t,s,s')z s ,ds'J 2 _, (6) 

Jo 



p(t, s, s') satisfy: 



dt 



g(t,s) 



d_ 

dt 



p(t,s,s') 



iuf(t,s)+2G(t)f(t,s) 

2iP(t,s) (7) 

iLug(t,s)-2F(t)f(t, S ) 

2F(t)g(t,s)+4G(t)f(t,s) 
2G{t)g(t,s)-2iP(t,s) (8) 

2iuip(t, s, s') + 2F(t)p(t, s, s') 



+ 2f(t,s)P(t,s , )-2g(t,s)P(t,s') (9) 
where the time-dependent functions F(t) = 
Jq a(t,s)f(t,s)ds, G(t) = f£ a(t, s)g(t,s)ds and 

P(t, s') = f Q a(t, s)p(t, s, s')ds together with a set 
of boundary conditions for /, g and p: f(s,s) = 1, 
p(t, s, t) = -ig(t, s), g(s, s) = 0, p(s, s, s') = 0. 

With the explicit O-operator |6]), the exact linear non- 
Markovian QSD equation may be compactly written into 
a time-local form, 

j t Mz) = ( - twJ, + J-Z t - F(t)J+J- - G{t)J+J z J„ 



i j P{t,s r )z s ,dsJ+J^Mz)- 



(10) 



Note that the time dependent coefficients F(t), G(t) and 
P(t, s') could be calculated once the correlation func- 
tion a(t,s) is explicitly given. It should be emphasized 
that the non-Markovian properties are encoded in a fi- 
nite width correlation function a(t,s), and so in the 
time-dependent coefficients F(t), G(t) and P(t,s') ap- 
pearing in the QSD equation (jlOl) . Clearly, the terms 
containing G(t) and P(t, s') give rise to the most im- 
portant correction to the Markovian dynamics. This 
can be easily seen if we take the Markov limit in which 
F(t) = 1/2, G(t) = P(t,s') = recovering the standard 
Markovian QSD equation. Eq. (jlOl) is the fundamental 
equation for the three-level system coupled to a bosonic 
environment described by @. 

For numerical simulations, it is more efficient to use 
the nonlinear non-Markovian QSD equation [l8[ for the 



where the coefficient functions f(t,s), g(t,s) and normalized state ipt 



obtained from Eq. (TTD")) : 



^-Jt = -uuJ z i> t + A t (J_)^ t + {J + ) t F{t)A t {J_)i> t - F{t)A t {J+ J_)j> t + (J+) t G(t)A t (J z J^ t 
at 

-G{t)A t {J+J z JJ)i) t + i{J + ) t [ P{t,s')z s ,ds'At(Jl)i>t-i [ P{t,s')z s ,ds'A t {J + J 2 _)j> t . (11) 

Jo Jo 



Here z t — z t + J* a* (t, s)(J + ) s ds is the shifted complex The correlation function a(t,s) at zero tempera- 



Gaussian process. 



ture is given by a(t,s) = £) k \g k \ 2 e 



2 p —iUk(t—s) 
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J °° dujS(w)e~ lu ( t ~ s \ where S(w) is the spectral density 
of the environment modes. The Eq. is capable of de- 
scribing non-Markovian dynamics with an arbitrary cor- 
relation function a(t, s), but for the sake of simplicity, 
we now consider correlation function a(t, s) — 5^e~ 7 ' t_s ' 
for the Ornstein-Uhlenbeck process. The Ornstein- 
Uhlenbeck process is a useful approach to modeling noisy 
relaxation with a finite environmental memory time scale 
I/7. When 7 — > 00, the environment memory time ap- 
proaches to zero and a(t,s) reduces to S(t — s), which 
corresponds to the Markov limit. Then Eq. (fTTj) reduces 
to the standard nonlinear QSD equation [J]: 



(a) 5 trajectories 



d j 



-iojJ z ^t + A t (J^){z t + (J + ) t )^ t --A t (J + J^)iJj t . 

(12) 

For the Ornstein-Uhlenbeck correlation, the partial dif- 
ferential equations Q, ([8]) and ([9]) could be converted 
into a set of ordinary differential equations: 



dF(t) 
dt 

dG(t) 
dt 

dP(t) 
dt 



7 



+ (-7 + iu)F(t) 



+2F(t)G{t) - 2£P(t), 
(- 1 + iuj)G(t)-2F 2 (t) 
+6F(t)G(t) - 2G 2 {t) - 2iP(t), 
-ilG(t)+2(-j + iuj)P(t) 
+AF{t)P{t) - 2G(t)P(t) 



(13) 



(14) 



(15) 



where P(t) — J Q a(t, s')P(t, s')ds' and the initial condi- 
tions are given by F(0) = G(0) = P(0) = 0. It is easy to 
show that P(t,s') = -iG(s')e&-r+ 2iu + 4F ^- 2G ^ ds . 
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FIG. 1. Ensemble average of (J) over 1000 realizations (red 
dot-dashed line for (Jx), green dashed line for (J y ) and blue 
solid line for (J z )) with different 7s. Here we choose uj — 1 
and the initial state \ip ) = 1/V~3(\0) + |1) + |2)). 

With the nonlinear non-Markovian QSD equation (jTTJ) , 
the simulations of the three-level system (a spin-1 par- 
ticle or a three-level atom) dynamics can be efficiently 
implemented by realizing Gaussian sample paths. We 
first calculated the ensemble average of the angular mo- 
mentum time evolution with different memory parameter 
7. The plots are shown in Fig. 1. The zero-temperature 
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(b) 1000 trajectories 




FIG. 2. Ensemble average of Purity = Tr[p 2 (t)] over different 
number of realizations, with red dot-dashed line for 7 = 0.2, 
green dashed line for 7 = 0.8 and blue solid line for 7 = 2.0. 
Here we choose uj = 1 and the initial state \ipo) — l/\/3(|0} + 
|1> + |2». 



environment prohibits the transitions from a lower level 
to a higher level, so for the three-level atom with an 
arbitrary initial state, the spontaneous emission always 
causes the system to decay into its ground level |0) in 
the long time limit. Therefore, when t — > 00, we get 
(Jx) = (Jy) — > and (J z ) — > —1, this scenario can be 
easily seen from the Markov limit shown in Fig. 1(d). 
The non-Markovian features of the environment for dif- 
ferent 7 are illustrated in Figs. 1(a), 1(b) and 1(c), where 
we can see that the transition of dynamics from non- 
Markovian to Markov regimes is dictated by environment 
memory time r = 1/7. Clearly, the non-Markovian fea- 
tures are lost when the system approaches to its Markov 
limit. The most important non-Markovian corrections 
are dominated by G(t) and the noise terms contained in 
the O-operator of Eq. ([6]), they become more significant 
in the case of smaller 7 (or longer r). As a consequence 
of long environmental memory, the dynamics of (J z ) has 
a long "tail", which means that it needs more time to 
reach its final steady state compared with the case of the 
Markov limit (Fig. 1(d)). The Markov dynamics emerges 
when r becomes shorter and shorter, that is, for the large 
7 ^> 1, G(t) and the noise terms can be effectively ne- 
glected. 

As a measure of degree of decoherence, we now con- 
sider the purity dynamics of the three-level system as 
shown in Fig. 2. It can be shown that the purity varies 
from 1 for a pure state to 1/d (d is the dimension of 
the density matrix, here we have d = 3) for a maxi- 
mally mixed state. We begin with a pure initial state 
\ipo) = 1/V3(|0) + |1) + 1 2)), so for the zero-temperature 
case, the final state of the three-level system is also pure. 
As shown in Fig. 2, the decoherence pathways of the 
three-level system are profoundly modified by the envi- 
ronmental memory parameter 7. An interesting feature 
of the quantum trajectory is that it can reveal how a 
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quantum state evolves into decoherence with only a small 
number of realizations. For example, Fig. 2(a) shows the 
results generated by only 5 realizations. Clearly, one can- 
not expect that the 5-realization simulations can quan- 
titatively reproduce an accurate description about de- 
coherence dynamics, nevertheless, they still give rise to 
a rather interesting qualitative picture about the gen- 
eral dissipative behaviors of the purity including the cor- 
rect information about the decay and recovery tenden- 
cies. For the numerical simulations with 1000 realizations 
shown in Fig. 2(b), it can be shown that a reliable phys- 
ical picture can be obtained. Specifically, when 7 = 0.2, 
the system exhibits a stronger non-Markovian feature, as 
such, the exchanges of quantum information (the energy 
distribution among the three levels |0),|1),|2) and the 
coherence between any two levels) with the environment 
via dissipation proceed slowly and information dissipated 
into environment may come back to the system in a finite 
time. Consequently, the decoherence time is effectively 
prolonged. On the contrary, for short memory time with 
7 = 2, as shown in Fig. 2(b), the system quickly evolves 
into a mixed state and then relaxes itself to the final 
pure state. When 7 = 0.8, the purity dynamics shows a 
moderate non-Markovian behavior. 

In conclusion, we present an exact non-Markovian 
quantum state diffusion equation for the dissipative 
three-level model described by Eq. @. We instigated 
significant and important progress by obtaining an ex- 
plicit form of the O-operator at zero temperature. With 
the time-local non-Markovian QSD equation, we are able 
to attack the transient property of quantum decoherence 
dynamics of the three-level system in all possible non- 
Markovian regimes and the well-known Markov limit. 
There are many important connections of our current 
work with multiple qubit and qutrit systems, a more de- 
tailed research into this subject will be useful. 
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